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An improvised technique of quintic 
hermite splines to discretize generalized 
Burgers—Huxley type equations 


I. Kaur, S. Arora*and I. Bala 


Abstract 


A mathematical collocation solution for generalized Burgers—Huxley 
and generalized Burgers—Fisher equations has been monitored using the 
weighted residual method with Hermite splines. In the space direction, 
quintic Hermite splines are introduced, while the time direction is dis- 
cretized using a finite difference approach. The technique is determined 
to be unconditionally stable, with order (h* + At) convergence. The tech- 
nique’s efficacy is tested using nonlinear partial differential equations. Two 
problems of the generalized Burgers—Huxley and Burgers—Fisher equations 
have been solved using a finite difference scheme as well as the quin- 
tic Hermite collocation method (FDQHCM) with varying impacts. The 
FDQHCM computer codes are written in MATLAB without transforming 
the nonlinear term to a linear term. The numerical findings are reported 
in weighted norms and in discrete form. To assess the technique’s appli- 
cability, numerical and exact values are compared, and a reasonably good 
agreement is recognized between the two. 
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1 Introduction 


Nonlinear partial differential equations nowadays turn out to be basic math- 
ematical methodologies to study multifarious structures like turbulence in 
fluid dynamics, convection-diffusion, flow through a shock wave roaming in 
viscid fluid, number theory, continuous stochastic processes, and so on. The 
diversity of physical phenomena in basic and applied sciences such as physics, 
chemistry, biology, computer science, electronics, and so on can be preem- 
inently described by these nonlinear equations. Various authors have used 
these nonlinear partial differential equations in different fields such as hy- 
drodynamics, solid mechanics, and so on; see [7, 9, 10, 11, 18, 25]. The 
Navier-Stokes equation is a fundamental fluid dynamics equation that may 
be simplified into a number of mathematical phenomena. By omitting the 
pressure factor, the generalized Burgers—Huxley equation (GBHE) simplifies 
this intricate equation. This equation explains how reaction, diffusion, and 
convection processes interact. Special instances of the GBHE include gen- 
eralized Burgers—Fisher equation (GBFE) and generalized Burgers equation. 
These equations are widely studied in fluid dynamics, gas dynamics, traveling 
wave solutions, and traffic flow. The one-dimensional GBHE, which refers to 
nerve pulse transmission in nerve fibers and waves in fluid crystals, can be 
written in the following form: 


Ou Pu 5 Ou 
u 


= the Q T\, 1 

Fp ep ge Hteoe + fu), (wt) € OK (0,7) (1) 
The initial and boundary conditions are presented as 

u(x,0) = g(a), (2) 

u(0, t) = Fit); and u(1, t) = fa(t), (3) 


where g(x), fi(t), and f2(t) are continuous functions in x and ft, respectively. 
Moreover, f(u) is a smooth function with a nonlinear nature that is defined 
on 2 x (0,7). In the theory of traveling wave solutions, it is significant. For 
f(u) = 0, equation (1) reduced to modified Burgers equation. On the other 
hand, in [24], it has studied the generalized Burgers equation with f(u) 4 0, 
although in a different way than Burgers-Huxley and Burgers—Fisher. For the 
GBHE, f(u) = Bu(1—u®)(u®—7) and for the GBFE, f(u) = Bu(1—u®), where 
Lt, 6, 8, and ¥ are real parameters. The viscosity factor ¢ is another distinction 
between the Burgers equation and Burgers—Huxley equation. The viscosity 
factor is typically assumed to be 1 in the Burgers—Huxley equation. However, 
it plays a vital role in turbulence theory for the Burgers equation. A specific 
nonlinear evolution equation is represented by the change in parameters. 
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For example, equation (1) presents the modified Burgers equation to explain 
wave propagation in nonlinear dissipative systems and various other physical 
contexts, such as sound waves in a viscous medium, where 3 = 0. Therefore, 
a significant study of GBHE with different case studies helps to analyze the 
behavior of different nonlinear equations from a wide perspective. 

Due to the wide applications of GBHE and GBFE, these equations are studied 
extensively by many investigators. Hammad and El-Azab [13] have used a 
collocation method with a 2N-order compact finite difference scheme, and 
Kushner and Matviychuk [17] proposed finite-dimensional dynamics to find 
the exact solution of Burgers—Huxley equation. Celik [8] proposed the Haar 
wavelet method, whereas Alharbi and Fahmy [1] have proposed an ADM- 
pade method to study the behavior of Burgers—Huxley equation. Javidi [16] 
followed a spectral collocation method to solve the GBHE, whereas Saha 
Ray and Gupta [21] followed a wavelet collocation for the same. The implicit 
exponential finite difference method has been followed in [14] to find the 
numerical solution to Burgers—Huxley equation. 

To solve the GBHE and GBFE, the present work proposes a quintic Hermite 
collocation with the forward finite difference technique (FDQHCM). Quintic 
Hermite collocation is a weighted residual approach with a Hermite basis, 
whereas finite difference scheme is the variational method. Both approaches 
are coupled to get numerical outputs that are compatible with numerical 
codes and are stable. 


2 Quintic Hermite collocation 


Quintic Hermite collocation method [4, 6, 5] is one of the weighted residual 
methods. Instead of Jacobi orthogonal polynomials, quintic Hermite inter- 
polating polynomials are used as the foundation function in this approach. 
Over the area 2, which is considered to be the union of intervals [x;_1, xj], 
the solution function is approximated by an approximation function that uses 
quintic Hermite interpolating polynomials as its basis. 

Hermite interpolating polynomials, which are an extension of Lagrangian 
interpolating polynomials, are of order 2k + 1, where k is a positive inte- 
ger. Hermite interpolating polynomials, on the other hand, are superior 
to Lagrangian interpolating polynomials in terms of applicability since they 
interpolate both the function and its kth-order derivative. Furthermore, La- 
grangian interpolating polynomials need a requirement of continuity at node 
locations that are not required by Hermite interpolating polynomials. As a 
result, quintic Hermite interpolating polynomials with k = 2 yield quintic 
Hermite interpolating polynomials that interpolate the function as well as 
its first- and second-order derivatives. Arora and Kaur [5] discussed the be- 
havior and structure of quintic Hermite polynomials in depth. Let u7(z,t) 
be the approximating function to be adjusted to eq. (1). Although the first- 
and second-order derivatives are interpolated by quintic Hermite interpolat- 
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ing polynomials, the boundary constraints are satisfied at boundary points. 
The approximating function is 


2 


u(x,t) = >> (Pi(w)ai(t) + Pi(a)bi(t) + Pi(a)ei(t)) (4) 


i=l 


where a;,b;, and cis are continuous functions of 't’ and P;, P;, and P; are 
smooth functions of ‘7’ and expressed as 


5 4 3 
g ( tite 2 15 ( Zoe +10 ( ene : BS eX Bea, 
Lji+1—Vj Lj+1—-Vj Ljp1—-Lj 


5 3 
P;(x) = 16 ( @-Lj—1 ) 15 ( uae ) 10 ( pa tyet ) 5. ga 4 <a< Ly, 


i-XLj-1 Di-LjZ-1 Li-XLj-1 
0 elsewhere, 
(wj41-2)° (wj41—2)* (aj41-2)° ; 
GiaesF iamee (ena? : a ee ee oe 
D. = L-Lj_ (a—aj~ @—2j-1) 
Pz) SG noe) Gea! Gan SPS ey 
0 elsewhere, 
(wj41-2)° (xj41—2)* (aj41-2)° : ; 
_ , as Gain . ora Lj SUS X41, 
D. = @—2;~-1) a LL; 
(2) US (rae FG carey e es ren Ti STS %;, 
0 elsewhere, 
where 
Pi(x3) = 6,8, 3 Ate = 0, £; (xj) =0, 45.9 = 1,2} 2056; 
Pi(zj)=0, P(a;)=4,, P, (vj) =0, i,7 =1,2,...,6, 


(x3) 

P,(x;) =0, (x5) = 90, P, (25) = 0%, ty 9 = Ly Qyaon5 6s 
The points xri,8 are the node points or mesh points, which are taken 
equidistant to get uniform mesh grid. The principle of collocation is ap- 
plied between two consecutive node points, that is, on [x;,2;+1]. The details 
are given in [2, 4, 5]. Therefore, to apply the principle of collocation, a new 
variable € is introduced within each sub-interval [x; x;+41] in such a way that 
as x varies from x; to 741, € varies from 0 to 1, and fh is the length of 
sub-interval [x; x;+41]. To reduce the complexity of system of equations, the 
interval length h is taken to be uniform. Therefore, after rearranging P;, P;, 
and P,, quintic Hermite polynomials take the form as given in Table 1. 


Table 1: Presentation of quintic Hermite splines 


i H; Hi Hi 

1 1—10€3 + 15é4 — 6€° —30€? + 60€3 — 30€4 —60E + 180€? — 120€? 

2 h(é—6€2 + 8é4— 36°) (1 — 1862 + 3262 — 15é4) (36 + 96E? — BOE?) 
3 Ace? 3634 3¢4 €5) BEng — 9g? + 108 bet) 4 (2 — 18¢ + 366? — 2069) 
4  -(e8 _ 264 + 5) h (3€? — 83 + 5é4) W (6€ — 24€? + 20€%) 

5 (10€3 — 15é4 + 6€5) (30€? — 60€3 + 30€*) (60€ — 180€? + 120€?) 

6 


h(—4€3 + 7€4 — 3€°) h(—12€? + 28€3 — 15é4) h(—24€ + 84€? — 60E3) 
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The values of quintic Hermite splines at x = 0 and x =1 are given in 
Tables 2 and 3. It is analyzed from these tables that these polynomial are 
either 0 or 1 at these points, which helps the approximating function to satisfy 
the Dirichlet or Neumann type boundary conditions. These quintic Hermite 
splines have the properties Hy(€) = Hs(1—&), H2(€) = —He(1 — §), and 
H3(€) = Ha(1 — €). 


Table 2: Presentation of quintic Hermite splines and its corresponding first- and second- 
order derivatives at €=0 


H, 1 Hi 0 HF O 
H, 0 Hi h HY 0 
H; 0 Hi, 0 HY Pp 
H, 0 Hi, 0 HY 0 
H; 0 Hi 0 HY 0 
He 0 Hi 0 HY 0 


Table 3: Presentation of quintic Hermite splines and its corresponding first- and second- 
order derivatives at €=1 


H, 0 Hi 0 HF O 
Hp 0 HS 0 HY 0 
H; 0 Hi 0 HY 0 
H, 0 Hi 0 HY Pp 
H, 1 HL 0 HY 0 
He 0 Hi h HY 0 


Therefore, equation (4) can be rewritten as 


6 
u7(é,t) = S> Hm(€)a%,@), (5) 


where H,,,(€)’s are quintic Hermite interpolating polynomials and a7, (t)’s are 
continuous functions of 't’ with 7 being the number of sub-divisions. To apply 
orthogonal collocation, zeros of fourth-order shifted Legendre polynomials 
have been taken as collocation points. These polynomials give better results 
on center as well as on average [2]. Therefore, the zeros of these polynomials 
have been chosen as collocation points. The diagrammatic representation of 
these splines is given in Figure 1. In Figure 2, the application of quintic 
Hermite splines has been shown. Further details of collocation points and 
the technique of quintic Hermite collocation are given in [5, 2]. 
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Figure 1: Diagrammatic behavior of quintic Hermite polynomials 


re Fh 6 a” og 


(y-1)™ edement (y+1)" dement 


(y)* cement 


Figure 2: Diagrammatic representation of quintic Hermite collocation scheme 


3 Implementation of FDOQHCM 


To implement the scheme of FDQHCM, equation (1) is discretized in time 
direction using forward finite difference scheme with step-size At 


Uj+1 — Us 


At = Uraj — [Us Un} se fis (6) 


where f; is f(w;). Initial and boundary conditions can be discretized in the 
following way at t = tj, uo(z) = g(x); u,;(0) = fi; and u;(1) = fo;, where 
fi; = fi(t;) and fo; = fe(t;). By simplifying, equation (6) converts into 


Uj+1 = Attire; _— Atuusur; + Atf; + U;. (7) 
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After applying quintic Hermite collocation on the variable u; in the space 
direction, equation (6) takes the following form: 


6 6 6 5 
Yana Ge) = “ H" (€)a7,(t3) — a Ss Hn(€)a(t)) 
m=1 ; m=1 

«(3 Hin(@alts)) +o Halts) + Dts, 


At the kth collocation point, equation (8) can be written as 


6 At 6 At 6 6 6 
De Hem Ona =F Dy Beans — =( » Hamas xD Akan, j 
m=1 m=1 


m=1 m=1 


6 
+ S- Hem, j + At fej; (9) 


m=1 


m=1,2,...,6, =2,3,4,5, 7 =1,2,...,m4, 


where H;,, is the mth interpolating polynomial at kth collocation point and 
Brm and Apm are, respectively, the second-order and the first-order deriva- 
tives of mth interpolating polynomial at kth collocation point. Also, Orn j 
is the mth collocation function at jth time step in 7” sub-domain and fr,j 
is the discretized function f(u) at the jth time step and the kth collocation 


point. 


Hence, after implementation of quintic Hermite collocation with forward 
finite difference scheme, collocation equations reduce to the following matrix 
form: 


Haj41 = Mya; + (Ha;)°(Moa;)+ f+ K, (10) 


where H = [Hm (fx)], Mi = [Bm (Ex) + Hm(Ex)], and Mz = [5°An(&)]. 
Moreover, K is obtained from the boundary conditions and f = [Atfy,,]. 
Also, G41 and G@; being the unknown collocation vectors at the (j + 1)th 
and jth time steps, respectively. In collective form, left-hand side coefficient 
matrix can be written as 


J = Mya; + (Haj)? (Moa;) + f+ K 
and A = [41]. 


Hence, in the combined form, equation (10) can be written as 


A=n-*y. (11) 
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Details of these matrices are given in [5, 3]. In the case of the first and last 
elements, the bandwidth is 4 x 5, whereas for the remaining elements, band- 
width is 4 x 6. Therefore, after combining all the collocation equations, a 
set of 4ne x 4ne equations appear, with ne being the total number of sub- 
divisions. The discretized set of equations has been solved in MATLAB. For 
initial approximation of Go, the initial condition of u at t = 0 has been taken 
at different collocation points. Then the loop was introduced to calculate a, 
Gg, and so on upto desired accuracy. The details of algorithm are mentioned 
hereunder. 

Algorithm: 

The algorithm of the solution technique is 

(i) Define the problem. 

(ii) Apply the finite difference technique in the time direction. 

(iii) Discretize the problem by applying quintic Hermite interpolating poly- 
nomials in space direction. 

(iv) Apply collocation points on the interpolating equations. 

(v) Solve the collocation equations using the following code. 

for t@=1:1:ti/j 

P(:, ii) = B— mux A.* (U%!!2) +U + betax j.*U. * (1— (U4). « ((Uset2) — 
gamma); 

Z(:, it) = inv(A) * P(:, ii); 

U0 = Z(:, it); 

end 


4 Stability analysis 


To study the stability of any partial differential equation, it is necessary that 
it should be linear. In the case of nonlinear partial differential equations, it 
is first linearized, and then the stability behavior is analyzed. Therefore, to 
study the stability of equation (1), it is first linearized by taking 0 = max(u®) 
and f(u) = Buf, where f = maa((u® —1)(y—u°)), for all (x,t) € Ox (0,T). 
As given in [20], a method is said to be A-stable if the region of absolute 
stability includes the region of Re(AAt) < 0. 

The region of absolute stability for FDQHCM is the set of all complex num- 
bers AAt such that 


Ju; | < C as t; > oo. (12) 
Let Aj, 7 = 1,2,...,ne, be the set of eigenvalues of the coefficient matrix 
a}, ;- The scheme of FDQHCM is said to be stable if Re(Aj, At) < 0 for all 


values of 1, 3,y, and 6. Then 
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6 i he 
De Am (Gini = FE ey Ain (Ging — HGP | Hin (Gin, 9 
m=1 m=1 m=1 
6 6 (13) 
+) Hal€)e gt Ae >, Halal f 
m=1 m=1 


The stability of the given technique can be checked on a linear equation. 
Therefore, equation (13) represents the quasi-linearized form of equation (8) 
to check the stability of the given technique. From Figure 3, it is observed 
that Re(A;, At) € [—1,0] for all values of y, 6,7, and 6, which justifies the 
stability of the numerical scheme. 


Imginary part of Eigen Values 
i=} 
1 


-1 1 1 1 
-1 -0.5 0 05 1 
Real Part of Eigen Values 


Figure 3: Behavior of eigenvalues for coefficient matrix given in equation (13) 


5 Convergence analysis 


The operator £ defined by = in spatial and time domains, is positive def- 
inite in L2(0,1), the space of all real valued Lebesgue measurable functions 
square integrable (0,1), for all t > 0. The definition given by [19] is quoted 
here: 

Consider a family of mathematical problems parametrized by singular pertur- 
bation parameter €, where € lies in the semi open intervalO0 <é¢ <1. Assume 
that each problem in the family has the unique solution denoted by ue and that 
each ue is approximated by a sequence of numerical solutions (Uz, Oe 4) 
where U,- is defined on the ON, representing the set of points in R and N 
is the discretization parameter. Then the numerical solutions U- are said 
to converge to the exact solution uz, if their exist a positive integer No and 
positive numbers C' and p, where No, C, and p are all independent of N and 
€, such that for all N > No, we have 


supo<e<i||U- > Ue| lan <CN?. 
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Here p is the rate of convergence and C’ is the error constant. 


Theorem 1 (Maximum Principle). [19] Let u(x) be a solution of an ad- 
vection diffusion equation with u(0) > 0 and u(1) > 0. Then L(u(z)) > 0 
for all x on the domain Q,L be the operator. Hence u(x) > 0 for all x in 9 


Lemma 1. [5, 3] Let H be the space of all Hermite interpolating polynomials 
of order 5 defined on the interval 0 < « < 1. Then 


6 
S- \Hm(x)| <6, — for all x € [0,1]. 


m=1 


Theorem 2. [12] Let H be the space of all quintic Hermite interpolating 
polynomials H of function u(x) defined on [a,b]. Then the rate of convergence 
of quintic Hermite interpolation on [a,b] is of order 6. Moreover, 


ju” (x) — H™ (x)|| < Cyn", =n =0,1,...,5 (14) 


where, the values of 7, are given in [12]. 


Theorem 3. [5, 3] Let U(z,t,;) be the quintic Hermite spline interpolate of 
u(a,t) from the space H of all Hermite interpolating polynomials of order 5 
defined on the interval 0 < x < 1 be the quintic Hermite spline interpolate 
of u(x,t), such that P,(x,t;) € C®({a,b]). Then the uniform error estimate 
is given by 

\|u(a, tj) — Ua, ty)Iloo < C(At + hi). (15) 


6 Numerical implementation 


To check the applicability of FDQHCM, numerical results have been com- 
pared to the analytic results. Stability of the proposed technique has been 
checked by L2-norm and L,,-norm as follows: 


L2-norm,  ||u||2 = |= >> wi (ul, t) — u7(&, t))?, (16) 


where wis represents the weight function corresponding to the collocation 
points and 


Loornorm, |}ullac = amax u(t) — u(Eis Oh (17) 


where u(&;,t) is the analytic solution and u7(&;,t) is the numerical solution 
obtained from FDQHCM. 
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Problem 1. Consider the GBHE given in equation (1) with f(u) = Gu(1 — 
u°)(u® — 7). The exact solution to the given equation is 


= 2 1/6 
u(x,t) = (2 +  tanh(ay(2 - a2t))) ; (18) 


where a] = Tarsy (— Hd +6/u? + 48(1+6)) and 


_ 2n7—(1+d—-F)(—Ht+V/ vw? +48(1+8)) 
oo 2145) : 
The initial and boundary conditions can be derived from the exact solution. 
From Tables 4—8, numerical values obtained from the FDQHCM have been 
compared to the exact values and of the results obtained from literature 


(13, 15, 23]. 


Table 4: Comparison of results for w= 1, 8 = 1, 6=1, and y = 0.001 


a ot Exact FDQHCM = Absolute Error [13] [23] [15] 

0.1 0.05 0.000500019 0.0005000110 7.9750x10-° 7.7006x10~° 7.72768x10~° 1.9372x 10-7 
0.1 0.000500025 0.0005000137 1.1300x107-8 =1.1268x 1078 1.12968x107-8 3.8743x 1077 
1 0.000500137 0.0005001205 1.6478x10~® 1.6863x10~-§ 1.68647x10—§ 3.8750x 1076 

0.5 0.05 0.000500069 0.0005000514 1.7607x10~8 = 1.7284x10~® 1.73534x1078 1.9373x 10-7 
0.1 0.000500075 0.0005000462 2.8837x10~8 2.8738x10~® 2.88305x10~8 3.8746x 10~7 
1 0.000500187 0.0005001405 4.6463x10~-°8 4.6841x10~® 4.68491x10—8 3.8753x 107 

0.9 0.05 0.000500119 0.0005001110 7.9760x10- 7.7006x10~° 7.72823x10—-9 1.9375x 1077 
0.1 0.000500125 0.0005001137 1.1301x10-8 1.1268x 1078 1.12980x10-® 3.8749x 10-7 
1 _0.000500237 0.0005002205 1.6481x10~® 1.6863x107~§ 1.68669x10~® 3.8756 x 1076 


Table 5: Comparison of results for 1 = 1, 8 = 1, and 7 = 0.001 


6=2 6=3 
2 ¢ Exact FDQHCM Absolute Error Exact FDQHCM Absolute Error 
0.1 0.05 0.0223614813 0.0223611210 3.60349x10—7 0.0793728109 0.0793714961 1.31479x10-§ 
0.1 0.0223617974 0.0223612705 5.26908x10—" 0.0793740199 0.0793720972 1.92270x10~® 
1 0.0223674857 0.0223666991 7.86609x10—" 0.0793957760 0.0793929065 2.86954x 107° 
0.5 0.05 0.0223634233 0.0223626133 8.09958x10—7 0.0793790070 0.0793760514 2.95560x10~® 
0.1 0.0223637393 0.0223623937 1.34559x10-® 0.0793802158 0.0793753057 4.91013x10~-® 
1 0.0223694271 0.0223672410 2.18610x10-® 0.0794019686 0.0793939936 7.97497x10~® 
0.9 0.05 0.0223653650 0.0223650047 3.60299x10—7 0.0793852021 0.0793838875 1.31464x10~® 
0.1 0.0223656810 0.0223651541 5.26903x10—" 0.0793864108 0.0793844882 1.92259x10~® 
1 0.0223713683 0.0223705817 7.86614x10—" 0.0794081601 0.0794052905 2.86962x10~°® 


It is observed that results obtained from the FDQHCM are at par with 
[13] and are better than [15, 23]. It is also observed that the FDQHCM gives 
error varying from 10~° to 10~° for varying values of py, 6,5, and 7. The 
results obtained from the FDQHCM appear to be more consistent than of 
[13]. In Table 9, L2-norm and L.-norm have been calculated for 4 = 6 = 1, 
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Table 6: Comparison of results for 1 = 1, 8 = 1, and 7 = 0.001 


6=4 


6=8 


2 ft Exact FDQHCM Absolute Error 


Exact 


FDQHCM Absolute Error 


0.1 0.05 0.1495399548 0.149537428 2.52679x 10~® 
0.1 0.1495423528 0.149538658 3.69481x10—§ 
1 0.1495854975 0.149579985 5.51247x10-® 
0.5 0.05 0.1495506669 0.149544987 5.67995x10—§ 
0.1 0.1495530645 0.149543629 9.43545x 1078 
1 0.1495961998 0.149580879 1.53208x 1075 
0.9 0.05 0.1495613768 0.149558851 2.52579x 10~® 
0.1 0.1495637738 0.149560080 3.69379x10—® 
1 0.1496068999 0.149601387 5.51287x 107 


0.38670979 0.386702959 6.83447x 10~® 
0.38671673 0.386706734 9.99399x 10-6 
0.38684140 0.386826506 1.48943x 1075 
0.38673162 0.386716259 1.53640x 10-5 
0.38673855 0.386713032 2.55228x 1075 
0.38686318 0.386821783 4.13948x 1075 
0.38675344 0.386746612 6.83197x10~® 
0.38676037 0.386750381 9.99200x 10~® 
0.38688495 0.386870055 1.489155x107° 


Table 7: Comparison of results for » = 1, 8 = 1, 6 = 16, and ¥ = 0.001 


x t Exact FDQHCM = Absolute Error 
0.1 0.05 0.6218689 0.6218574 1.14208x10~° 
0.1 0.6218811 0.6218644 1.66988x10-° 
1 0.6221000 0.6220752 2.48186x10-° 
0.5 0.05 0.6218956 0.6218699  2.56702x10~° 
0.1 0.6219078 0.6218651 4.26434x107° 
1 0.6221265  0.6220576 6.89814x10~° 
0.9 0.05 0.6219223 0.6219108  1.14113x107° 
0.1 0.6219344 0.6219178  1.66897x10-° 
1 0.6221531 0.6221283 2.48109x10-° 


Table 8: Comparison of results for w= 0, 8 = 1, 6=1, and 7 = 0.001 


ot Exact FDQHCM = Absolute Error 


[13] 


[23] [15] 


0.1 0.05 0.000500030 0.00050001989 
0.1 0.000500043 0.00050002762 
1 0.000500268 0.00050024509 
0.5 0.05 0.000500101 0.00050007775 
0.1 0.000500113 0.00050007495 
1 0.000500338 0.00050027582 
0.9 0.05 0.000500172 0.00050016131 
0.1 0.000500184 0.00050016904 
1 0.000500409 0.00050038651 


1.0286 x10— 

1.5044 10-8 
2.2466 x 10-8 
2.3132x10-8 
3.8429x 10-8 
6.2443 x 1078 
1.0283 x 10-8 
1.5047 10-8 
2.2464 10-8 


1.0269x 10-8 
1.5027 10-8 
2.2488 x 10-8 
2.3049 10-8 
3.8324 1078 
6.2465 x 1078 
1.0269 10-8 
1.5027 10-8 
2.2488x 10-8 


1.0303x10—® 1.8747x 10-7 
1.5063x10—8 3.7493x 107-7 
2.2488x10-8 3.7500x 10° 
2.3138x10—8 1.8749x 10-7 
3.8441x10—8 3.7498 10-7 
6.2465 10-8 3.7504 10-8 
1.0303x10—8 1.8751x 1077 
1.5063x 10-8 3.7502x10—7 
2.2488x10-8 3.7509x 10~—® 


¥ = 0.001, and for 6 = 1,2,3, respectively. Both the norms are lying within 


the range of 10-8 to 107°. 


Table 9: La-norm and Loo-norm for p= 1, 8 = 1, and 7 = 0.001 


t Lz -— norm 


Doo-norm 


6=1 6=2 


6=3 


6=1 


6=2 6=3 


0.1 6.328299x10—9 2.954021x10~7 1.077846x10-® 2.880650x10-® 1 
0.25 9.224090x10-° 4.305708x10—7 1.571002x10~-® 4.271075x10-® 1.993759x10-® 7.274619x10~® 
0.5 1.000496x10~® 4.669971x10—7 1.703806x10~® 4.645850x10~® 2.168583x10~® 7.911991x10~° 


.344790x 10-© 4.906733 x 10-5 
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From Figures 4—7, three-dimensional behavior of u(x,t) has been pre- 
sented for u = 8 = 1, 7 = 0.001, and 6 = 1,2,3,4, respectively, which shows 
that the values are lying between 0 and 1 for time ranging from 0 to 0.01. 


Figure 5: Three-dimensional behavior of u(x,t) for » = 1, 8 = 1, y = 0.001, and 6 = 2. 


Problem 2. Consider the GBFE given in equation (1) for ¢ = 1 with f(u) = 
Bu(1—u°). The exact solution for the GBFE is 


1 1 1/6 
u(x,t) = (5 + 5 tanh(a3 (a — ast))) ; (19) 
where a3 = wy and a4 = see des The initial condition and boundary 


conditions can be obtained from the exact solution. From Tables 10-17, 
numerical values obtained from the FDQHCM have been compared to the 
exact values and the results obtained from [13, 23]. It is observed that results 
obtained from the FDQHCM are at par with the results obtained from [13] 
and are better than [22]. However, in Tables 10 and 17, results by [13] are 
slightly better than the FDQHCM, but this is overcome by the simplicity of 
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Figure 7: Three-dimensional behavior of u(x,t) for » = 1, 8 =1, y = 0.001, and 6 = 4. 


the FDQHCM as compared to the technique of former. The absolute error 
varies from 10~!! to 107° for varying values of pz, 3, and 6. In Table 18, L2- 
norm and L,,.-norm have been calculated for uw = 6 = 0.001 and 6 = 1,2,3, 


respectively. 

Table 10: Comparison of results for = 0.001, 6 = 0.001, and 6 = 1 

et Exact FDQHCM Absolute Error [13] [22] [15] 

0.1 0.001 0.499988 0.499988 9.999x10~!9 5.8147x10~1! 1.01 x10~7 9.68763x 10° 
0.005 0.499989 0.499989  6.999x10-9 2.6058x10~19 4.38 x10~" 1.93753x10~§ 
0.01 0.499990 0.499990 1.100x10~8 4.4599x1071° 7.53 x1077 1.93752x107° 

0.5 0.001 0.499938 0.499938 4.460x10~1% 5.6241x10711 1.04 x10~7 9.68691x10~® 
0.005 0.499939 0.499939 9.310x10~!3 3.0621x10~19 5.21 x10-" 1.93738x10~§ 
0.01 0.499940 0.499940 9.985x10-19 6.1867x10~1° 1.04 x10-® 1.93738x 1075 

0.9 0.001 0.499888 0.499888 9.980x10-1° 5.8135x10-1!! 1.01 x10-7 9.68619x10~® 
0.005 0.499889 0.499889  6.998x10~-9 2.6053x1071° 4.38 x1077 1.93724x107§ 
0.01 0.499890 0.499890  1.100x10~8 4.4591x1071° 7.53 x10~7 1.93724x 1075 
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11: Comparison of results for 4 = 0.001, 8 = 0.001 


6=2 


6=3 


x t Exact FDQHCM 


Absolute Error 


Exact FDQHCM Absolute Error 


0.1 0.001 0.70709535 
0.005 0.70709676 
0.01 0.70709853 

0.5 0.001 0.70704821 
0.005 0.70704962 
0.01 0.70705139 

0.9 0.001 0.70700106 
0.005 0.70700248 
0.01 0.70700424 


0.70709535 
0.70709677 
0.70709855 
0.70704821 
0.70704962 
0.70705139 
0.70700106 
0.70700249 
0.70700426 


4.65367 x 10 19 
1.02299x 1078 
1.54395 1078 
2.5057 10719 
7.9174x10711 
1.2627x 10714 
1.17716x10~9 
1.07532x1078 
1.57272x10-8 


0.79369100 0.79369100 5.3640x10~1° 
0.79369259 0.79369260 1.1100x10-§ 
0.79369457 0.79369459 1.7813x10~§ 
0.79365131 0.79365131 4.9914x10~1° 
0.79365290 0.79365290 9.5943x10~14 
0.79365489 0.79365489 4.1912x10~1° 
0.79361162 0.79361162 1.4294x10~9 
0.79361321 0.79361322 1.1675x10~§ 
0.79361519 0.79361521 1.6992x10~§ 


Table 12: Comparison of results for » = 0.001, 6 = 0.001, and 6 = 4 


x t Exact FDQHCM = Absolute Error [22] 
0.1 0.001 0.84088843  0.84088843  1.37598x10-9 1.75 x10-8 
0.005 0.84089011 0.84089012 1.25399x10-8 7.37 x1077 
0.01  0.84089221 0.84089223  1.82591x10-® 1.27 x10~° 
0.5 0.001 0.84085479 0.84085479 2.0918x10-1!9 1.75 x1078 
0.005 0.84085647 0.84085647 = 1.7133x107-! 8.77 x 10-7 
0.01 0.84085857 0.84085857 6.3821x107!° 1.75 x107® 
0.9 0.001 0.84082114 0.84082115 1.0779x 1079 1.75 x10-8 
0.005 0.84082283 0.84082284 1.1838 1078 7.38 x10-7 
0.01 0.84082493 0.84082495 1.8053 1078 1.27 x107§ 
Table 13: Comparison of results for 4, = 0.001 and 6 = 0.001 
6=8 6=16 
x t Exact FDQHCM Absolute Error Exact FDQHCM Absolute Error 
0.1 0.001 0.91699941 0.91699941 8.421940x10-!9 0.9576009 0.9576009  1.023x10~9 
0.005 0.91700124 0.91697903 2.222702x10-5 0.9576028 0.9576029 4.177x10~-8 
0.01 0.91700353 0.91695864 4.491036x10—° 0.9576052 0.9576053 4.778x 1078 
0.5 0.001 0.91697903 0.91700124 2.221422x10-5 0.9575897 0.9575897  2.435x10-8 
0.005 0.91698086 0.91698086 3.297860x 10-19 0.9575916 0.9575916 2.357x10~—1° 
0.01 0.91698315 0.91696048 2.267573x10—° 0.9575940 0.9575940 1.435x1078 
0.9 0.001 0.91695864 0.91700353 4.488876x10—5 0.9575784 0.9575784 6.336x10-° 
0.005 0.91696048 0.91698315 2.266323x10—5 0.9578803 0.9575803  2.227x10-8 
0.01 0.91696277 0.91696277 1.947571x10~® 0.9575827 0.9575827 2.092x10~8 
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6=2 
[13] & [22] 


6=4 


FDQHCM FDQHCM [13] 


0.1 0.1 1.9313x10-5 1.76638x10—5 1.121x10~-® 5.9667x10~© 1.26230x10-5 1.343x10—5 
0.3 7.0917x10—° 2.51379 x107° 1.600x10-> 2.1779x10-° 1.79797 x107° 1.919x107 
0.5 1.2556x107° 2.61751 x107-° 1.667x10-5 3.8540x107° 1.87212x107° 2.001x10~ 

0.5 0.1 6.4389x10~® 4.49179 x10—5 2.904x10-5 2.0158x10~® 3.21358 x10—5 3.489x 10-5 
0.3 5.0445x107-°5 6.91014 x107° 4.468x10—5 1.5604x107° 4.94694x107-° 5.373x107° 


0.5 1.0433x10-4 7.24595x10—5 4.687x10~5 3.2214x10-5 5.18702 x10—5 5.641x10~5 
0.9 0.1 1.9529x10-5 1.75391 x10-5 1.154x10-5 6.1067x10~® 1.25646x10-5 1.393x10~5 
1.7920 x107-5 1.981x1075 
0.5 1.2709x10-4 2.60482 x10-5 1.711x10~5 3.9456x10-° 1.86610 x10—5 2.065x10—5 


0.3 7.1767x10—5 2.50111 x10 1.643x107° 2.2305x107° 


Table 15: Comparison of absolute errors for w = 0.1, 8 = —0.0025, and 6 = 8 


a ¢t FDQHCM [13] [22] 

0.1 O01 4.8295x10-" 7.65875x10-® ~=1.471x10-5 
0.3 1.5645x107® 1.09129x107-5 2.107x107 
0.5 2.6947x10-® 1.13630 x10-5 = 2.203« 1075 

0.5 0.1 1.8448x107" 1.95143 x1075 = 3.832x 1075 
0.3 1.1484x10-® 3.00445 x10-5 5.911«1075 
0.5 2.2738K107® 3.15019 x107> 6.218x1075 

0.9 0.1 4.9302x10-’ 7.63553 x1075 = 1.533x1075 
0.3 1.6078x10-® 1.08887 x107-5 2.183« 1075 
0.5 2.7772x107® 1.13383 x107> 2.280107 


Table 16: Comparison of absolute errors for » = 1 and 6 = 0 


6=3 56=8 
x t FDQHCM [13] FDQHCM [13] 

0.1 0.0001 2.3559x10-™ 5.7870x10-"! -8.59409x10-"! —-2.7729x10-" 
0.0005 2.6056x10-® 2.0870x10—-5 1.35615x107®8  1.0782x10-5 
0.001 6.7400x10-" 4.5521x1075 3.86425x1077  2.3520x107° 

0.5 0.0001 4.6311x107!9 5.1671x1071!! 9.37x1071! 2.0409 x 1071! 
0.0005 2.6853x107!® = 1.7895x10-5 =4.82443x10-!9 9.4066 x 10~® 
0.001 4.8165x107-!9 4.0262x1075 2.49371x107!9 =. 2.1163x 107° 

0.9 0.0001 4.5141«107-!° 5.3143«107!!  1.40942x10719 = 1.8090x 10714 
0.0005 2.9347x1078 =1.6546x1075 ~—-1.7873x107-8 ~—-8.7305x 1078 
0.001 7.5936x10-" 3.6150x10—-5 4.92825x10-"  =1.9075x 1075 
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Table 17: Comparison of absolute errors for » = 1 and 6 = 1 


6=1 6=2 

x t FDQHCM FDQHCM [13] [22] 

0.1 0.0001 1.2969x10-© —3.03026x10-! ~—4.4788x10-™ 1.55 x10-> 
0.0005 5.9286x10-5 9.30236x10-8 2.4700 x10-5 7.62 x10~5 
0.001 1.9592x10-4 = 1.62366x 107° ~=—5.3872 x107° 1.50 x 1074 

0.5 0.0001 5.9153x107!® 1.99555x1071!9 7.1278 x10-!!—- 1.83 «1075 
0.0005 1.5519x107-!® 2.98812x10-!° 2.0953 x10-5 9.14 x1075 
0.001 9.1902x10~-8  2.91233x107!® 4.7151 «1075 ~=—-1.83 «1074 

0.9 0.0001 1.1642x10-® 2.34912x107!® 5.5348 x107!9 2.07 x107® 
0.0005 5.3457x10-° =: 9.77826x107-® ~=1.9314 x107-5 ~—-:1.02 x 10-4 
0.001 1.7683x10-4 = 1.70078x107® 4.2211 x10~5 2.00 x10~4 

Table 18: Lo-norm and L.o-norm for ~ = 0.001 and 6 = 0.001 

t Lz — norm Doo-norm 

6=1 6=2 6=3 6=1 6=2 6=3 


0.01 6.996123x10~° 8.760110x10~° 9.773331x10~° 4.937500x10~® 6.633549x10—5 7.056646 x 10—§ 
0.05 1.243733x10-® 1.776084x10—8 2.003356x10—* 6.687500x10-® 1.007198x10—7 1.096476x10-7 
0.1 1.680553x10~® 2.385368x10—8 2.735586x10—® 7.375000x10—® 1.090978x10—7 1.229183x10~7 


Both the norms are lying within the range of 10~° to 10~". From Figures 
8-11, three-dimensional behavior of u(x,t) has been presented for uw = 6 = 
0.001 and 6 = 1,2,3,4, which shows that the values are lying between 0 and 
1 for time ranging from 0 to 0.01. 


Figure 8: Three-dimensional behavior of u(a,t) for u = 0.001, 8 = 0.001, y = 0.001, and 


6=1. 
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Figure 9: Three-dimensional behavior of u(a,t) for 4 = 0.001, 8 = 0.001, y = 0.001, and 
6 =2. 


u(x,t) 


Figure 10: Three-dimensional behavior of u(z,t) for 4 = 0.001, 8 = 0.001, y = 0.001, 
and 6 = 3. 


Figure 11: Three-dimensional behavior of u(z,t) for 4 = 0.001, 8 = 0.001, y = 0.001, 
and 6 = 4. 
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7 Conclusions 


The technique of the FDQHCM has been implemented successfully on the 
GBHE and GBFE. The technique is a combination of the weighted residual 
method and the finite difference method, which gives stability to the numeri- 
cal results and is easily adaptable to computer codes. Numerical results have 
been calculated for a vast range of parameters, even for those parameters not 
been calculated earlier. The technique is found to be more consistent than 


[13, 


15, 22, 23]. It was also shown that the FDQHCM could be applied to 


nonlinear partial differential equations of higher order too. 
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